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Abstract 



Doppler tracking of interplanetary spacecraft provides the only method 
presently available for broad-band searches of low frequency gravitational 
waves (~ 10~^ — IHz). The instruments have a peak sensitivity around 
the reciprocal of the round-trip light-time T (~ 10^ — 10^ sec) of the radio 
link connecting the Earth to the space-probe and therefore are particularly 
suitable to search for coalescing binaries containing massive black holes in 
galactic nuclei. A number of Doppler experiments - the most recent involving 
the probes ULYSSES, GALILEO and MARS OBSERVER - have been car- 
ried out so far; moreover, in 2002-2004 the CASSINI spacecraft will perform 
three 40 days data acquisition runs with expected sensitivity about twenty 
times better than that achieved so far. 

Central aims of this paper are: (i) to explore, as a function of the relevant 
instrumental and astrophysical parameters, the Doppler output produced by 
in-spiral signals - sinusoids of increasing frequency and amplitude (the so- 
called chirp); (ii) to identify the most important parameter regions where to 
concentrate intense and dedicated data analysis; (iii) to analyze the all-sky 
and all-frequency sensitivity of the CASSINI's experiments, with particular 
emphasis on possible astrophysical targets, such as our Galactic Centre and 
the Virgo Cluster. 

We consider first an ideal situation in which the spectrum of the noise 
is white and there are no cutoffs in the instrumental band; we can define 
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an ideal signal-to-noise ratio (SNR) which depends in a simple way on the 
fundamental parameters of the source - chirp mass Ai and luminosity distance 
- and the experiment - round-trip light-time and noise spectral level. For 
any real experiment we define the sensitivity function T as the degradation 
of the SNR with respect to its ideal value due to a coloured spectrum, the 
experiment finite duration Ti, the accessible frequency band {fb, fe) of the 
signal and the source's location in the sky. We show that the actual value of 
T crucially depends on the overlap of the band (/b, /e) with the instrument 
response: the sensitivity is best when l/T and fe coincides with the 
frequency corresponding to the beginning of the merging phase. Furthermore, 
for any ft, and Ti, there is an optimal value of the chirp mass - the critical 
chirp mass Aic oc /^^ '^^^ - that produces the largest sensitivity function; 

lower values of M. correspond to a smaller bandwidth and less SNR. Also the 
optimal source's location in the sky strongly depends on (/;,, fe). 

We show that the largest distance at which a source is detectable with 
CASSINI's experiments is ~ 600 Mpc and is attained for massive black holes 
of comparable masses ~ 10"^ Mq and fb ~ 10"^ Hz. Sources not far from coa- 
lescence in the Virgo cluster with 10® Mq ^ Ai ^ 10^ Mq would be detectable 
with a SNR '-^ 1 — 30. The SNR and the range of accessible masses reduce 
drastically when a smaller mass ratio is considered. We then turn the atten- 
tion to galactic observations, in particular on the detectability of a coalescing 
binary in the Galactic Centre, where a small black hole of mass M2 could be 
orbiting around the central massive one Mi ~ 2 x 10^ Mq. CASSINI would 
be able to pick up such systems with M2 > 50 Mq; for M2 ^ 10^ Mq the SNR 
could be as high as ~ 100 — 1000. It may also be possible to detect such bi- 
naries in more than one of the three CASSINI experiments, thus re-enforcing 
the confidence of detection. 



PACS numbers: 04.80.Nn, 97.60.Lf, 98.62.Js 
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I. INTRODUCTION 



The search for gravitational waves has greatly grown in the last few years: four ground 
based laser interferometers - LIGO 0, VIRGO 0, GEO600 g and TAMA [| - with optimal 
sensitivity in the band ~ 10 Hz — 1 kHz are now under construction and are scheduled to 
be in operation by the turn of the century; at the same time, the sensitivity of acoustic 
narrow-band devices in the kHz regime is steadily improving [^j-^]. In this high frequency 
band the most interesting sources are connected to collapsed objects with mass in the range 
~ 1 — 1000 Mq. However, at low frequencies (below a few Hz) detectors on Earth are severely 
impaired by the seismic noise. To observe sources of higher mass one must use instruments 
in space, sensitive in the mHz band. 

The Doppler tracking of interplanetary spacecraft is the only method presently 
available to search for gravitational waves in the low frequency regime (~ 10~^ — IHz). 
Several experiments have been carried out so far, all of them with non dedicated space- 
probes: VOYAGER 1 in|, PIONEER 10 and 11 [0,0], ULYSSES GALILEO and 

They suffer from the limitations of instrumentation designed for 
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other purposes and from several disturbances, including those due to the ground sector; 
however, their sensitivity is astrophysically interesting, even though the rate of detectable 
events is uncertain and probably low, given the current instrument performances. At present 
a large amount of data is available: a major experiment lasted 28 days was carried out with 
ULYSSES in February/March 1992 and a 20 days coincidence experiment was performed 
in the Spring 1993 with the spacecraft GALILEO, MARS-OBSERVER and ULYSSES 
(see Table 1). 

The space-probe CASSINI - a NASA/ESA/ASI joint mission [l7| - represents the next 
step in Doppler experiments: launched on October 15th 1997 with primary target the study 
of the Saturn system, the spacecraft carries on board much improved instrumentation and 
will perform three long (40 days each) data acquisition runs in 2002, 2003 and 2004 to search 
for gravitational waves |]18| with expected sensitivity about twenty times better than that 
achieved so far. The next far away and much more ambitious projects will probably involve 
interferometers in space with arms of millions of km [ p!U|pU[ (see also |2T[] and references 
therein for a review on detectors in space). 

In the low frequency band we expect five main types of sources (for an extensive discussion 
regarding the whole gravitational wave spectrum we refer the reader to [p2| , p3[ ): 



Catastrophic, wide-band collapses of large masses, possibly in dense concentrations of 
matter and stars, leading to the formation of a massive black hole (MBH); the search 
for such signals motivated the original proposal of Doppler experiments [pi]. 



2. Short-period binary systems of solar mass compact objects [pB|-p7[; among them the 



binary pulsar PSR 1913+16 [p8| , p9| , whose study over the past 25 years has provided 
the most clear indirect evidence of the existence of gravitational waves ; however, 
stellar mass binaries are not detectable with present and near-future Doppler experi- 
ments, due to inadequate sensitivity and short duration. 



3. Solar mass compact objects orbiting a massive black hole in galactic cores [31-33 
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4. Binary systems of massive black holes spiraling together toward their final coalescence 

MM- 



5. Backgrounds of gravitational waves of primordial origin or generated by the superpo- 



sition of radiation coming from unresolved binaries ||37|-p|. 



Here we will consider signals emitted by binary systems containing massive black holes 
(MBH's), therefore sources at point 3 and 4, and the main goals of the paper are the explo- 
ration of the structure of the Doppler output produced by in-spiral signals from coalescing 
binaries and the analysis of the all-sky/all-frequency sensitivity of Doppler experiments with 
emphasis on the CASSINI mission. 

Compelling arguments suggest the presence of MBH's in the nuclei of most galaxies, as 



the near-inevitable byproduct of the infall of gas in their potential well |l42|-^ and MBH's 
are invoked to explain a number of phenomena, such as the activity of quasars and active 
galactic nuclei p5H48[]. However, the observational evidences of their existence come mainly 
from observations of relatively nearby galaxies, whose nuclei do not show significant activity. 
Central dark masses - that are believed to be "dead quasars" |^ - are mainly inferred by 
studying the spatial distribution and velocity of gas and/or stars in galactic cores. About 
ten candidate massive black holes, with mass in the range ~ 10^ — 10^ Mq, are known today 



and we refer the reader to ||50|-^ for recent reviews. The most striking evidences come from 
the mapping of the gas motion via the 1.35 cm water maser emission line in NCG 4258 and 
the observations in the near infrared band of the star motion in Our Galactic Centre. In the 
case of NCG 4258, observations with the Very Long Baseline Array have allowed to measure 
the velocity of the gas in the disk with an accuracy of 1 km/sec: the disc rotates with a 
velocity exactly in agreement with Kepler low, due to the presence of an obscure object of 
mass ~ 3.6 x 10^ |Q. Regarding Our Galactic Centre, recent observations with the 
New Technology Telescope of the proper motion of stars in the core have shown that their 
speed scales as the inverse of the square root of the distance from the center up to 20000 
km/sec; the inferred mass of the dark body is ^ 2.5 x 10^ Mq |55]. In both cases only exotic 



(and highly implausible) alternatives to a MBH can still be considered p6| , p7 

Many galaxies have experienced at least a merger since the epoch z > 2 [^,^]. If a 
MBH is present in each core of both interacting galaxies, the two massive objects would fall 
into the common potential well, forming a pair which losses energy and angular momentum 
through dynamical friction; eventually a binary forms and driven by radiation reaction 
precedes toward the final coalescence. This scenario provides, therefore, a quite natural 
way of producing massive binary black holes (MBHB's) [^|. The apparent precession and 



bending of jets in active galaxies P^j60|| and Doppler shifted broad emission line peaks 
in quasars |6l[] are possible indirect evidences of the existence of such systems. Moreover 
observations of MBHB's in the nuclei 1928+738 |2[, OJ287 ||,|6|] and 3C 390.3 |5[ have 
been claimed. 

MBH's in the center of galaxies are surrounded by a dense cluster of ordinary stars, whose 
tidal disruption provides the gas to refuel the central object, white dwarfs, neutron stars 
and low mass back holes, probably in the mass range ~ 5 — 100 Mq. Unless ordinary stars, 
compact objects are not disrupted by the hole and can occasionally be captured, presumably 
on highly eccentric orbits, forming binary systems with life shorter than the Hubble time. 
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Moreover the orbits of objects with mass ^ SMq around the central massive body would 
not be perturbed by other stars in the core providing therefore clean gravitational wave 
signals, that can be used to extract valuable information about the central object p6| , |67[ . 



Taking the data analysis point of view, where the structure of the signal determines 
the choice of the optimal algorithm, it is convenient to divide the entire process of binary 
coalescence in three separate phases: 

(a) Adiabatic in-spiral: the bodies are driven by radiation reaction from their initial 
distance to smaller separations, losing energy and angular momentum. The emitted gravi- 
tational wave corresponds to a sinusoid of increasing frequency and amplitude (the so-called 
chirp). In this paper we will make two important assumptions: (i) when the signal enters 
the sensitivity window of the detector, the orbit has already been circularized; in fact, the 
eccentricity e diminishes with the frequency / according to e^(t) oc f~^^^^{t) |6^. A binary 



born with a period of 100 years, say, even with e ~ 1 would carry a residual eccentricity 
^ 10~^ when the signal can be picked up by the instrument. For binaries composed of 
two MBH's, that have undergone a common evolution inside a galactic core, this assump- 
tion is reasonable, even-though the full understanding of massive binary evolution before 



radiation reaction takes over still lacks [3S-71]. The condition e <^ 1 is likely to be vio- 



lated for solar mass compact objects orbiting a massive one, where the secondary body is 
captured in a highly elliptic orbit, which maintains a not negligible eccentricity throughout 
the in-spiral [0. This would bring about a different scenario, with shorter life-time and 



emission of gravitational waves at multiples of the orbital frequency ||7^. (ii) In the mod 



elling of the waveform we will use the lowest-order Newtonian quadrupolar approximation, 
as we are mainly aiming at the estimation of the signal-to-noise ratio (SNR) and not at the 



construction of signal templates for matching filters [^3|,^ . 

(b) Merger, when the body separation r reaches a value of about 6 M, where M is the 
source's total mass, the orbit becomes unstable and the binary begins the final deadly phase 
of plunge-in. The signal consists in a wide-band burst of duration ~ M, but the details of 
the emitted radiation are in large part still unknown. 

(c) Ring-down: the single black hole formed during the merging phase settles down 
oscillating according to quasi-normal modes. The emitted waves consist of a superposition 
of exponentially damped sinusoids. 

In this paper we will concentrate on radiation emitted during the in-spiral phase. Such 
signals can be usefully classified according to the best algorithm to extract them in a given 
data set iTBI: 



I. Periodic signals: the rate of change of the frequency is too small to be detected with 
the instrumental frequency resolution. Here the ordinary techniques for the search of 
periodic oscillations of unknown frequency apply. 

II. Linearly "chirped" signals: the frequency drift is larger than the frequency resolution 
of the detector but the change in frequency due to its second time derivative during 
the observation time can be neglected, so that we can assume a linear increase of f{t) 
with a constant rate /3 = df/dt. An efficient technique to search for these signals 



has been developed in ||75|-[77[] and already applied to several data sets [0]; for every 
allowed /3, the chirping signal is reduced to a periodic one and standard methods can 
be applied. 
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These two classes correspond to sources which are far from coalescence - in fact df{t)/dt oc 
see Eq. Q - and, for a given distance, weaker; therefore, the search depth (the 
maximum distance at which a binary would be detectable) is very limited: in |]T5[ we showed 



that, with signals of this class, the Galactic Centre is marginally accessible with the best 
data available today. 

III. General chirps: the assumed signal has the correct evolution in frequency, but the 
merger does not occur inside the record. 

IV. Chirps and bursts: the record contains the final part of the chirping waveform, and 
the wide-band burst produced just before the engulfment of one black hole into the 
horizon of the other. 

These classes of signals can be distinguished in the record by the frequency and its change; 
for a given frequency, its drift increases with the class. Roughly speaking, this ordering 
corresponds to stronger sources and, therefore, to larger attainable distances. 

Although in the final phase of the in-spiral the two bodies have high velocities {v/c> 0.1) 
and the description of their gravitational wave emission requires a relativistic approach 
we shall use the quadrupole Newtonian approximation for the estimation of the signal-to- 
noise ratio (SNR). The use of this approximation has the advantage that the in-spiral signal 
depends on only one physical parameter, the chirp mass Ad, see Eq. (^^ ). In principle it 
can be derived from the measurement of the frequency / and its rate of change df/dt; then 
the amplitude A oc Ai^^^ / D provides the distance D of the source. We can therefore 
determine, for any given experiment, the largest attainable distance; we call it, for a given 
class, the ken of the experiment. This is another example, similar to the one pointed out by 
Schutz in W^ , of the power of gravitational wave astronomy when a simple source model is 



available. Central aim of the paper is to fully develop the richness of these concepts when 
applied to Doppler experiments and to link possible astrophysical sources (like those in the 
Virgo cluster and/or in the Galactic Centre) to the actual data. In particular we will apply 
the main results to the set of experiments scheduled during the CASSINI's mission. 

Doppler tracking of interplanetary spacecraft provides a wide band detector of gravita- 
tional waves in the low frequency regime (but it can also be used as narrow-band xylophone 
instrument |]T9[). A very stable electromagnetic signal in the GHz radio band is continuously 
transmitted from the Earth to the spacecraft and coherently transponded back in order to 
monitor the change of their relative velocity. The optimal sensitivity of the detector is at 
frequencies of the order of the inverse of the round-trip light-time T of the Earth-spacecraft 
radio link. However, since in practice the duration of the experiment is not large, only a 
limited frequency interval of a chirping signal is accessible. There are two modes of search 
for a signal: if the frequency interval is small, narrow band search, one looks for sources 
far from the final plunge-in, hence at a small distance. In a wide band search the signal is 
stronger and one can detect radiation emitted by binary systems further away. In the past, 
the analysis of Doppler data has been confined mainly to the narrow band case (classes I and 
II); in this paper we concentrate on the more interesting wide band case and aim at sources 
at larger distances, but detectable only for a smaller fraction of their lifetime. The largest 
"ken" is, of course, attained for class IV, which merges into the class of generic, wide-band 
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bursts. The search for bursts in the records of the 1993 GALILEO-MARS OBSERVER- 
ULYSSES coincidence experiment is currently in progress |16||. However this analysis must 



be a rough one, as the structure of the signal in this strongly relativistic regime is essentially 
unknown and the use of data collected in coincidence by several instruments is crucial. 

The the paper is organized as follows. In Sec. ^ we review the basic concepts and 
formula (partially to fix notation) regarding the gravitational wave emission from binaries. 
In Sec. fT| we introduce the signal y{t) produced at the output of a Doppler detector 
by a gravitational wave train h{t): in the Fourier space, y{f) is the product of h{f) and 
the Doppler three-pulse response function reif)', we review the main properties of rg{f) 
and in particular we stress its different structure in the low (/ ^ ^/T) and high (/ ^ ^/T) 
frequency regimes. Sec. contains the main results of the paper: a thorough analysis of 
the signal-to-noise ratio produced by in-spiral signals in Doppler experiments as a function 
of the parameters that characterize the source and the instrument: the chirp mass Ai (and 
the mass ratio), the instantaneous emission frequency of the signal when the instrument 
is "turned on", the duration of the observation Ti, the structure of the instrument noise 
spectral density 5'„(/), the location 9 of the source in the sky with respect to the detector 
arm and the round-trip light-time T. This analysis is carried out in full generality and can 
be applied to any Doppler experiment. We define the ideal SNR, pid, as that corresponding 
to unlimited bandwidth, 6* = 7r/2 and flat noise spectrum: depends only on A4 ,T, D and 
the noise spectral level 5*0. Of course, in a real experiment several factors - among which 
the finite observation time, the signal bandwidth, the frequency-dependent noise spectral 
density and the position of the source in the sky - contribute to degrade the SNR. Indeed, 
we define the sensitivity function T as the actual SNR with respect to its ideal value. T 
is extensively studied in Sec. [IV A| . In particular we show that the maximum value of T 
is achieved when 0.1/T ^ fi, ^ 1/T and the binary chirp mass is close to what we call its 
critical value Aic oc Ti ^^^ f^^^^ . In Sec. p.V B| we study the CASSINI's sensitivity (and 
compare it with that of ULYSSES), with emphasis on possible sources in the Virgo Cluster 
and the Galactic Center. We show that the maximum reachable distance is ~ 600 Mpc; 
binaries close to the final coalescence with 10^ Mq ^ M. ^ 10^ Mq and comparable masses 
would be detectable in the Virgo cluster with SNR up to 30; in galactic searches, CASSINI 
would be able to pick up (possibly in all three data sets) signals from the in-spiral of a 
secondary black hole of mass M2 ^ 50 onto the central one Mi ^ 2 x 10® Mq. Finally, 
Sec. V contains a discussion about the probability of success of Doppler experiments based 
on simple conventional astrophysical scenarios and our conclusions. 

We take units in which c = G = 1, so that velocities are dimensionless and 10® Mq ~ 
4.926 sec. 



II. THE EXPECTED SIGNAL 

A. The chirped precursor 

We consider a binary system of two compact objects of mass Mi and M2; M = Mi -|- M2 
and p = M1M2/M are the total and reduced mass, respectively; the orbital parameters 
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evolve secularly due to the loss of energy and angular momentum by emission of gravitational 
waves. 

We shall assume that, when the gravitational signal enters the sensitivity window of the 
detector, the radiation reaction has already circularized the orbit. Under this assumption, 
in the Newtonian quadrupolar approximation, the frequency / of the emitted gravitational 
wave is twice the orbital frequency of the binary, namely 



fit) 



M 



7rV3(t) 



(2.1) 



where r{t) is the orbital separation. In the Newtonian model, f{t) evolves according to ||80 



/(t) = ^^8/3_M5/3/n/3(i) 
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where the dot indicates the time derivative and 



iv 



M 



<l/4) 



(2.2) 



(2.3) 



is the so-called chirp mass; in this approximation Ai is the only dynamical parameter that 
regulates the evolution of the frequency and the amplitude of the wave. 

Let us consider a signal emitted by a binary at a luminosity distance D and let l be the 
angle between the wave propagation direction and the orbital angular momentum L. The 
strain at the detector reads 



h{t) = (t) cos 2ip + (t) sin 2ip 



A{t) cos 



271 f dt'f{t')+^o 



(2.4) 



where h+{t) and (t) are the two independent polarization states 0], and ip a polarization 
angle; 



A//5/3 

A{t) = 2Q{i,v)^-[nf{t)f' 



D 



1.3 X 10-'''Q{i,ip) 



' D ' 
Mpc 



/ 



2/3 



M 



5/3 



lO-^Hz 



(2.5) 



is the amplitude of the signal and 



(fQ = tan 



2 cos t sin(2(y9) 



COs(2(y9)(l + COS^ l) 

gives the polarization phase. The quantity 

Q'^{L,ip) = [cos^(2v9)(l + cos^t)^ + 4cos^isin2(2v9)] 



(2.6) 



(2.7) 



depends on the orientation of the binary and, for random values, (Q^) = 8/5, where () 
stands for the average with respect to l and (p: Q{l, ip) ranges from (for l = 7r/2, (f = 7r/4) 
to 2 (for 6 = or vr) 
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Integrating Eq. ( p.2|) , it is straightforward to derive the frequency evolution of the wave: 



1/5x8/8 / .X -3/8 



and therefore the phase of the gravitational signal reads 



t -t\^/^ 



m = 2nl rft7(0 = <^>n - 2 (^^j ; (2.9) 
equivalently, 

t(/)=tn-5(87r/)-8/3Al-5/3^ (2.10) 



$(/) = $„- 2(87rA^/)-^/^ (2.11) 

From Eq. ( |2.2| ), one can also derive the number of Newtonian wave cycles spent by a binary 
while it sweeps the relevant frequency interval 

mf)= ldf = Af^--{8nMfr'/\ (2.12) 

Jf f TT 

tn, $n and Afn are integration constants, defined as the values that t{f), $(/) and A/'(/) 
formally take when / = 00 (r = 0). Of course, the signal must be cut off when the in-spiral 
phase ends and the merger begins. The fully relativistic two-body problem is still unsolved; 
it is mainly investigated using post-Newtonian approximations (see [|8T| , |82[| and references 
therein) and numerical techniques (see P^-pB[ and references therein). For simplicity we 
will neglect the radiation coming from frequencies higher than 

= 63757m - " " (10^) " • ("^' 

corresponding to the innermost stable circular orbit risco = 6 M for a test mass in the 
Schwarzschild field of a mass M. For r < risco (and velocity larger than fisco — 0.408) orbits 
are unstable. For general black holes and mass ratios we can write risco = Qk M, with /c ?s 1 
(several methods have been used to estimate the exact value of k, but a satisfactory answer 
still lacks; see p7| - |100|| for further details). In the Newtonian approximation we assume that 
/isco = f(tc) separates, at the time tc, the in-spiral phase from the broad band plunge-in. 
The major feature of the relativistic corrections is, therefore, to anticipate the final plunge. 
If an experiment takes place in the time interval from t = to t = Ti, we define 

ft^fit = 0) (2.14) 

as the frequency at which the signal enters the record at t = 0; the frequency at the end, 
after the observation time Ti, is 



fe = f{t = Ti) = min 



Jb \ t / ' 



(2.15) 
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If — Ti/t„)~^/^ > /isco the record includes the merger signal and is of class IV; otherwise 
we have class III signals. The Newtonian time to coalescence for a system radiating at 
frequency at the beginning of the record is given by: 

r„ = 5(87r/,)-8/3^-5/3. (2.16) 

It is also useful to note that the time for a signal to sweep from the frequency fb to the 
beginning of the final coalescence, that we suppose to occur at / = /isco, is 

= 5 (8vr)-«/3A<-V3 (;-8/3 _ ^^-^8/3) _ (2.17) 

We define 

/ T \ '^^^ 

" y + 5(8vr/.eo)-'«/^ -M-^/^j ^'-'^^ 

as the initial frequency of a gravitational wave whose final frequency is /isco! it depends on 
M., 7] and Ti and corresponds to the largest initial frequency that produces class III signals. 

On the base of the model of the waveform we can now formalize the distinction of the 
signals in classes that we have introduced in Sec. I. Being based upon the actual record, 
it is particularly useful for data analysis. It refers to the frequency resolution 1/Ti of the 
experiment and will allow us to divide the two-dimensional parameter space {Ai, ft), where 
the in-spiral signal is defined, into a periodic region, a linear region and a general chirp 
region (of class III and IV, depending whether the final burst is within the record or not). 
This distinction is justified by the great simplicity and low computational load of the data 
analysis when the frequency is constant (class I) or varies linearly with time (class II); these 
two cases are extensively discussed, and applied to the data of ULYSSES' second opposition, 



in 151. When 



5 

we have periodic signals; when 



T2/(t) = ^8/3j.2_^5/3y.ll/3 < j (^j^^^ (2.19) 



^ 16896 ^16/3 ^3 ^10/3 fl9/3 < 1 J.2 y^^) > ^ (^j^^^ jj^ (2.20) 

2 23 



linear signals; 



>1 and A < /isco (class III/IV) (2.21) 



correspond to general chirps ||101|| . Referring to a chirp mass A4 and observation time Ti, 



the initial frequency fb that separates out the four regimes reads: 



fb 



5.1xlO-(Tjfe,)"^^"(To^)"^^"Hz (/-//) 
5-4 X 10-^ {j^y'^'' (lo^r'"^'' Hz (// - ///) . (2.2 



1.9 X 10-3 (4,)3/^ 1 + 8.7x10^ {4,r (icfc) (lo^) " 



-8/3 



Hz {III -IV) 
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The partition of the plane {J^, fh) in classes for a single CASSINI's experiment (Ti = 
40 days) is given in Fig. |l]. The values of r„ that correspond to the transition between two 
adjacent classes read 

2.6xl0^(^)-^/"(^)"/"sec (//-///) . (2.23) 
[l.lxl03(^)-^/^(^)sec {HI -IV) 



As we have mentioned in the Introduction, class I to IV correspond to gravitational wave 
emission increasingly closer to the final merger and, therefore, to stronger signals. In Fig. 
H we show the frequency and the characteristic time r„ that separate the four classes as 
function of the observation time. For typical present and future data sets (see Table |) and 
reasonable sources, most signals would show up in class III or IV. 

We introduce now the Fourier transform h{f) of the real function h{t): 

Kf)= e-^-'f'h{t)dt =h*{-f) (2.24) 



and hereafter we will work in the more convenient Fourier space. When the frequency f{t) 
and the amplitude A{t) (Eqs. ( p.2[ ) and ( p.5| )) vary over a time scale much longer than 
1/ fit) it is convenient to use the stationary phase approximation in order to compute the 
integral ( |2.24| ). In our case 



where we have neglected numerical coefficients of proportionality. The evaluation of the 
actual error involved in this approximation is a delicate matter; however, in the limit 

[-M/(t)]'^'«l, (2.26) 

we expect this approximation to be satisfactory for the assessment of the signal-to-noise 
ratio. Note also that 

[Mf{t)f' < (A</.co)'/' oc (^^ j oc , (2.27) 



so that the limit (|2.26| ) is formally fulfilled for a small mass ratio. Of course, as usual, the 



approximation will be fearlessly pushed to the limits of validity. 

In the stationary phase approximation to the integral (|2.24|) , for a frequency / only a 
small interval 5t contributes around the time t{f) where 

f{t)=f. (2.28) 

Note the slight inconsistency of notation: / is now the independent variable in the Fourier 
space, while earlier we have introduced f{t) as a function of time. The time interval is of 
order 
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(2.29) 



corresponding to a local bandwidth 6f = l/6t. A finite record in (0, Ti) will pick up only the 
part of the spectrum between fb and fe] since dropping the data outside (0, Ti) is equivalent 
to convolving with sin(/Ti) / /, one wonders how this can kill the Fourier components outside 
{fb, fe)- The answer is that outside the frequency range defined by / = fit), {0 < t < Ti), 
the phase of the Fourier integral (|2.24|) has fast oscillations. The Fourier amplitude is then 
of order 



ASt 



D 



f 



-7/6 



The full expression of the Fourier transform of h{t), Eq. 
stationary phase approximation is ||102||: 



(2.30) 

4]), computed according to the 



Hf) 




/</. 
/>/i. 



(2.31) 



where 



A 



is the amplitude and 



1/2 



47r2/3 



7W5/6 



^(/) = 2nftn - $n 



^ + li8nMf)-'/' 



(2.32) 



(2.33) 



the phase. 

The existence of an analytic expression of the in-spiral signal in the Fourier space justifies 
its use. As compared with the analysis in the time domain, it offers two main advantages: 
the instrumental response function, see Eqs. (|3.1| )- (|3.4| ), is a simple multiplicative factor 
and we can easily deal with coloured noise (see Sections |ITT| and [TVI ) . 



B. The final burst 

The theoretical knowledge of the signal emitted during the merger of two black holes is 
at present quite poor. It involves the computation of the solutions of the Einstein equation 
in highly non linear regimes and this issue is currently tackled by several groups mainly 
by means of intensive numerical computation on state-of-the art machines (see p4|-pm and 
references therein). Here we present a simple model for the energy spectral density dE/df 
emitted during the merger, following the approach presented in ||103|| ; although quite naive 



and crude, it is reasonable for an order of magnitude comparison of the sensitivity of Doppler 
experiments to in-spiral and merger signals, which shall be carried out in Sec. |V[ 

The wave emitted during the final stages of the life of a binary system is spread over 
a wide band. For sake of simplicity, and lack of more detailed knowledge, we assume it to 
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be confined between the final frequency of tlie in-spiral, i.e. /isco, and the characteristic 
frequency /qnm at which the ring-down emission sets in. We will assume 

/^m ^ = 3.2 X 10-2 ( f.. ] Hz (2.34) 

that roughly corresponds to the frequency of the quasi-normal oscillation for the I = m = 2 
mode of a Kerr black hole with dimensionless spin parameter ~ 0.98. This choice of /qnm 
can be justified as follows: the I = m = 2 mode, very likely to be excited at the end of the 
black hole plunge-in, is that with the longest damping time and the resulting massive black 
hole will quite probably carry a large intrinsic angular momentum. 

The total energy -Emerg radiated during the merger phase, and confined into the frequency 
interval (/isco, /qnm), is proportional to /i^/M ||104|| and can be written as 



£^mcrg = emcrgi^(j^) M, (2.35) 

where emerg denotes the efficiency of the process and F{fj,/M) = {Afi/MY; of course, a 
precise estimate of emerg would require a full numerical calculation, but a reasonable order of 
magnitude guess yields 

Cmerg ~ 0.1 (scc ||103| and references therein). Lacking any evidence 
about the structure of the energy spectrum, we will make the simple assumption that -Emerg 
is uniformly distributed over the band (/qnm, /isco), so that 

J rp ( -Emerg /■ <^ f <^ f 

"^merg ) (f^^^^f,^^^) iisco ^ J ^ Jqnm _ ^2.36) 



df [0 / < /isco and / > /, 



qnm 



We recall here that the energy spectrum of the radiation emitted during the in-spiral phase 
is given by 

dE^ _ / ^ M'l' r'/' f < /isco _ (2.37) 



^^/ I / > / 



ISCO 

isco 



Eqs. (|2.36| ) and ( 2.37 ) will enable us to compare the detectability of coalescing binaries 



during the in-spiral and merger phase with Doppler experiments, as the optimal signal-to- 



noise ratio is simply related to dE/df, see Eq. ( [4.5|) ; this analysis is carried out in Sec. 



IV 



III. DOPPLER RESPONSE TO IN-SPIRAL SIGNALS 

In a Doppler experiment, the Earth and a spacecraft are used as the end-points of a 



gravitational wave detector [|I0[. A radio link is transmitted from the Earth to the spacecraft, 
coherently transponded and sent back to the Earth, where, at the time t, its frequency is 
measured with great accuracy; comparing the emitted and received frequency - uq and i/(t), 
respectively - one determines the Doppler shift (i^(t) — i^o)/^o as a function of time. The 
contribution to this effect produced by a gravitational wave h(t) is usually called the Doppler 
signal s{t) and reads [M 
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s{t) = / dt're{t-t')h{t') ; (3.1) 



008 6* — 1 ,, cos6' + l \ l + cos6'^, , , 

re{t) = 6{t) - cos 5 ft t\ + 6{t - T) (3.2) 

is the characteristic three-pulse response of the instrument, as the incoming signal h{t) is 
repeated in the detector output at three different times (and with different amphtude): t, 
t — (cos 6* + l)T/2 and t — T; rg(t) depends on the angle 9 between the spacecraft and the 
source and the round-trip light-time T = 2L out to the distance L of the probe. In the 



frequency domain the Doppler signal (|3.1| ) takes a very simple form: 

~s{f) = Mf)Hf) , (3.3) 
where h{f) is given by Eq. (|2.31| ) and the Fourier transform of ro{t) reads 

feif) = ^ cos^exp + cose)fT] + ^^exp[2mfT] . (3.4) 

Since the frequency / appears multiplied by the round-trip light-time T, we introduce the 
dimensionless variable 

e^/T. (3.5) 

In agreement with the Weak Equivalence Principle, the response of the instrument goes 
to zero in the low frequency limit ^ ^ 1 (this is the regime in which ground and space-based 
interferometers work): 

rg{^) = iTc^sin'^ 9 - ^n^^^ sin'^ e{2 + cose) + 0{^^) . (3.6) 

To lowest order the angle 9 only affects the coefficient. As ^ increases, rei^) develops 
modulations in amplitude and, more importantly, in phase, that depend on 6. The beats 
produced by the superposition of the signals at the times t, t — (cos 6* + l)T/2 and t — T 
result in minima and maxima of the amplitude superimposed to the general increasing trend 
A(t) oc P^^{t), see Eq. ([^.5]), as shown in Fig ^; they strongly depend on the angle 9 and 
provide a sensitive signature for its determination. 

The transversal character of gravitational waves is reflected in the forward {\9\ <^ 1) and 
backward {\9'\ = Ivr — ^ 1) limits, when the wave travels in a direction almost parallel to 
the Earth-spacecraft line: 

reiO = J [(1 + 2vrzOe^'^^« - l] + 0(9') m « 1) , (3.7) 
f,(^) = ^ [e2-« - 1 + 2n^] + 0(9'^) i\9'\ = \n - 9\ <^ 1) ; (3.8) 
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both tend to zero; they have the same hmiting amphtude, but different phases. Note that 
the two hmits are not uniform and are violated at high frequency ^ 1). 
For generic values of 6 it is convenient to use the squared modulus of (p.4|): 

, ,n 3cos^^ + l cos^6' — 1 , ^, 
\re{0? = ^ + ^ cos(2vrO 

— cos ^^(cos 9 + 1) cos[(cos 9 — l)7r^] 

- cos ^(cos ^ - 1) cos[(cos 9 + 1)71^] , (3.9) 



which, in the limit <^ 1, reads 



l_±(l_±cos2^)7r2e2 



+ 0{e). (3.10) 



I^e(0p5 which is shown in Fig ^, is an even function of cos^, sum of four functions periodic 
in ^, with frequencies 0, 1,(1 — cos9)/2 and (1 + cos^)/2, respectively. Its average is 
(3cos^^^ + l)/2. As example, for 9 7r/2 the response reads 

|fe(Or = sin2(7rO + [(l-cos(7rO)' 

-27ii,m{'ni)]{9-^-f + 0[{9-'^t] (3.11) 

and it has a maximum for = 1/3. For large frequencies it develops more and more lobes 
with the angular scale 1/^ (see Fig. ^). Averaging |r0(^)p over the direction 9 one obtains 
the mean square response of the detector: 



\re{0\' = W \fe{0\'d{cos9) 



2J-1 



which, in the low and high frequency limit, reads: 



cos(27rO + — ^ sin(27rO , (3.12) 



f,-r{i-inr)+oie) (^«i) .3^3^ 

l-icos(27rO (e»l) ■ 



During an experiment neither the angle 9 nor the round-trip light-time T are exactly 
constant, mainly due to the motion of the Earth; their time- variation induces changes in the 
phase of the filter of order ttA9 fT and tt/AT. For the latter AT ^ 2 V(sTi/c ^2x IQ-^Ti 
[v^ ~ 30 km/sec is the orbital velocity of the Earth) and the effect is relevant when 

2xlO"V/Ti^l, (3.14) 

while the angle 9 affects the phase by a smaller amount. For CASSINI's experiments this oc- 
curs at about 5 mHz, which is above the most important sensitivity region of the instrument 
(see next Section). Indeed, we will consider constant both T and 9. 
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IV. SENSITIVITY OF DOPPLER EXPERIMENTS 



The Fourier output of the detector 

y{f) = ~sif)+n{f) (4.1) 

is the sum of the signal s(/) and the noise n{f)- The noise is assumed to be stationary, with 
zero mean and spectral power density 

{n*{f)n{n) = S{f-nSM)- (4.2) 

We consider single-sided noise spectra Sn{f), formally from to +oo, and take into ac- 
count the finite instrumental bandwidth with suitably large values at the boundaries of the 



sensitivity window. When Sn{f) is white the corresponding Allan deviation ||105 | 




(4.3) 

is inversely proportional to the square root of the integration time r; for coloured spectra, see, 
e.g., ||106|| . For simplicity we use a continuous spectrum 5'„(/), although, over an experiment 



of duration Ti and with sampling time Ag, the spectrum consists of discrete Fourier modes 
spaced by l/Ti, from 1/Ti to the Nyquist frequency l/(2As). 

Since the shape of the signal in the in-spiral phase can be accurately predicted, the tech- 
nique of matched filtering is particularly suitable (see, e.g., |]73| , [7^ and references therein). 



It can be shown that if the signal s{f) is known to within a constant (and real) amplitude 
h, its least square estimator has a variance given by 

A)^,..4rf^<«/. (4.4) 

The quantity p is appropriately called the signal-to-noise ratio and is related to the energy 
spectrum of the radiation dE/df by the relation 

2 1 r MfW dE 



It is useful to introduce an ideal reference value of the SNR and to discuss its dependence 
on the main parameters. From Eq. ( |3.13 ), we can approximate the instrumental filter rg, 
as follows: 

< ^) (4 6) 

(e>l) • ^^-^^ 

With a white noise spectrum Sn = cr^r = const., and taking into account the frequency 
dependence of the signal (see Eqs. ( p.31|) , ( p.3|) and ( [1.6|) ), we have 

|s(/)P f'-^ (^<i) ,,,, 

^ \ f-V3 ^ . (4.7) 




Snif) I ^ (e > 1) 

y 
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This is a decreasing function of ^, both sides are integrable, so that in this ideal case the 
bulk of the sensitivity comes from / ~ 1/T (cfr. also Fig |^). We can define the ideal 
SNR by taking ^f, = /f,T <^ 1, = /g T ^ 1 and the direction angle 6 = 7r/2, for which 
1^^71/2(0 P = sin^(7r^); from Eqs. (|2.31|) , (|2.32|) and (|4.4|) , averaging over the source angles 

L, ip, we get 

Pm = K———, (4.8) 



where 



K = 



The ideal ken - the largest distance at which a source is ideally detectable at pid = 1 - is 

(4,10) 



For the CASSINI's mission (see Table |), with T = 10^ sec and a^r = 9 x lO-^^Hz"^ , it 
has the value 



'id 



45 (Tri-) L.,.. Mpc. (4.11) 



VlO^Moy VlO%ec/ V9xlO"^^Hz 



For ULYSSES the the ideal ken for Ai = IO^Mq, is ~ 3 Mpc; in the coincidence experiment, 
the values for GALILEO and MARS-OBSERVER are ~ 0.5 Mpc and ~ 2 Mpc, respectively, 
where the values of cij^, r and T have been chosen according to Table |. 

The ideal value ([4.8|) for the SNR neglects several degradation factors, that are actually 
present in real experiments, including: 

1. The direction of the source may be unfavorable. 

2. The frequency interval {fb, fe) swept by the signal is finite because the record has a 
finite length and the instrument has its own cutoffs. 

3. The noise spectrum is generally red, which damages the SNR in the critical low fre- 
quency band. 

Only experiments can provide a safe estimate of the noise structure, especially when aim- 
ing at high sensitivities or wide frequency windows. Past Doppler data can be approximated 
with a simple analytical expression made up with power law contributions: 

5n(/) = 5oE(|) ^ , (4.12) 
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where 



So = S„(l/r) = al{T)T 



(4.13) 



sets the spectral level; the frequencies fj and the exponents aj, fulfilling 



3 



1, 



(4.14) 



characterize the different contributions to the noise. We can define the transition frequencies 
/i2 and as the intercept between two power approximations; for example, the transition 
frequency /12 between region 1 and 2 (characterized by the index ai and 02, respectively) 
is given by 



In general, high frequencies (say, above ~ 0.2 — 1 Hz) are dominated by thermal noise, with 
the index 03 = 2; in the case of ULYSSES, this spectral region is also perturbed by the 
spacecraft rotational dynamics, so that the data have been actually low-pass filtered to 5 x 
10~^ Hz = /23 (for ULYSSES we can therefore assume 03 1 and we have formally achieved 
that by setting 03 = 10). The frequency window between /12 and /23 can be characterized 
by a power law with index a2 ~ —0.5, due to propagation noise [p!^p!5| . The analysis of 
past experiments has been mostly confined to this band, in particular to frequencies larger 
than ^ 1/T = /12; in the case of ULYSSES, the procedure for the generation of residuals 
essentially cuts off all frequencies below ~ 10~^ Hz; we have formally achieved this by 
adopting ai = —10 (see Table |I|). 

As we have already anticipated in the Introduction (and we shall show more in detail 
in the next Sections), one of the main conclusions of this paper is that low frequencies (i.e. 
^ 1/T) are crucial to achieve large SNR and to widen the range of accessible masses. The 
noise spectrum for fT ^ 1 is to a large extent unknown, as systematic effects, mostly related 
to the orbit determination process and the tropospheric correction, are likely to come into 
play. For the goal of assessing the expected SNR for CASSINI we take the view that these 
systematic effects do not exceed the propagation noise above 10""^ Hz and keep 0:2 = —0.5. 
The level Sq is given by the main specification for the Allan deviation, i.e., cXy = 3 x 10~^^ at 
10"^ sec, corresponding to S',i(10~'*Hz) = 9 x 10~^^ sec; we also limit our band to a minimum 
frequency of 10~^ Hz. For the interval up to fu = 10^^ Hz the only experimental evidences 
come from ULYSSES data, which showed an approximate behavior as and, accordingly, 
we chose ai = —2; this value, however, is conservative, as in the processing of ULYSSES 
data no particular attention was taken at the reduction of the noise below 10~^ Hz. The 
effect of possible different values of ai will be explored in Sec. [IV A\ Thermal noise, with 
as = 2, can be expected to set in at about 0.1 Hz; however, the limited sampling time will 
prevent to go much beyond. In Fig. ^ we show the noise spectral density of CASSINI and 
ULYSSES, computed using Eqs. ([4.12|) and ( |4.13| ) with the parameters given in Table 

Keeping Sq and the frequencies fj (j=l,2,3) fixed, we now write the SNR in the form 




(4.15) 
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p' 



SoD^ 



(4.16) 



The sensitivity function 




-7/3 



^(1 



measures the change of SNR with respect to the "ideal" case, Eq. 

T(0,cx),7r/2;0) = 1. 



and fulfils 



(4.17) 



(4.18) 



Generally this function is less than unity and gives a degradation; we must however point 
out that the value 7r/2 for the angle 6 is not always (i.e. for any arbitrary choice of aj) the 
best; furthermore, we have defined the "white noise" considering the level of the real noise 
at / = 1/T, so that Sn{f) may actually be smaller than the white noise for / > 1/T. Notice 
also that when ai > —2/3 the SNR spectral density is not integrable and the total SNR is 
dominated by the low frequency end fb- 

The crucial physical parameter of a binary is its chirp mass Ai, which characterizes the 
rate of change of the gravitational wave frequency and sets the relationship (except for class 
IV) 



M 



(87r)8/5 



ry-i — 3/5 8/5 

-'i Jb 



'-il 



8/3 



3/5 



Mr 



1 



Je 



8/3' 



3/5 



(class I - III) . 



(4.19) 



We see that, given and Ti, there is a critical value of the chirp mass 

53/5 



(87r)8/5 
9.2 X 10^ 



rji — 3/5 r — 8/5 

^ 1 h 



-8/5 



-3/5 



10-4 Hz 



40 days 



(4.20) 



above which no measurement is possible, in class I-to-III. At this critical value, when -C 
/isco; ^1 basically coincides with the time to coalescence r^. 

When M. <^ M.C1 the signal is intrinsically narrow-band, and so are the measurements; 
the fractional bandwidth 



A 



/e ~ fb 

h 



is small and given by 



A 



3 /A^ _ 3 Ti 
8 \Mj ~ 8 7> 



(4.21) 



(4.22) 
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This case, in which there is httle power and information, is less interesting. 

When the data set includes the merger (class IV), instead of Eq. ( 4.19| ) we have 



M = M, 




'A 

, /isc 



8/3N 



3/5 



(class IV) . 



(4.23) 



The coefficient within square brackets can be greater than unity, so that the inequality 
< A^c does not apply. 

Since, with decreasing frequency, the spectral density of the SNR decreases or has a 
weaker growth below 1/T, it is very convenient to have T/^ < 1. When = feT -C 1 the 
response ivr^ sin^ Eq. ( |3.6| ), weakens the signal. When > 1 the response introduces 
strong modulations that depend on 6 and we have full sensitivity. This is the interesting 



region and will be explored in Sec. P^V A| using the dynamical behavior of the source. 

Before discussing in detail the behavior of the SNR in observations of in-spiral signals as 
function of the relevant parameters, it is instructive to compare the sensitivity of Doppler 
experiments to the radiation emitted during the in-spiral and merger phase. If pinsp and 
Pmerg corrcspoud to the optimal SNR achievable by observations of the same source in the 
two regimes, from Eq. ([4.5|) we can write: 



Pinsp 
P 



mcrg 



/is CO 
B 



df 



df 



/qnm \rg (/)|2 dEn.crM) 



/is, 



df 



df 



1/2 



(4.24) 



We stress that in the previous relationship we have assumed optimal filtering not only for 
the in-spiral but also for the merger signal; if for the former the assumption is realistic, cfr. 
| |107| |, for the latter it is indeed very optimistic. Notice also that pinsp is the largest SNR 
obtainable for a chirp, as we integrate over the widest frequency band (/s, /isco) for a data 
set of length Ti (cfr. Sec. |1V Aj) , and we have averaged the three-pulse response function 



over the angle 6, cfr. Eq ( |3.12| ). For a given source, if (pinsp/Pmerg) > 1 (< 1) it is therefore 
'more convenient' to search for it by means of observations of the in-spiral (merger) signal. 
Inserting Eqs. (g]3|) and into Eq. (142^ ), the ratio of SNR's becomes 



Pinsp 
Pmerg 



I.IJ 



0.1 



-1/2 



T 



10^ sec 



1/6 



M 



1/3 



-1/2 



2 X W Mq 



5 X 106 M©^ 



(4.25) 



where 



J 



'(T/b) 



3) Jfgg)l_ 



d^ 



J (T /isco) eSr^iO 



d^ 



1/2 



(4.26) 



measures the relative overlap of the in-spiral and merger signals with respect to the instru- 
ment response. In Fig. |^ we show (pinsp/ Pmcrg) and J as a function of the chirp mass for 
two different mass ratios, 1 and 0.01, for a nominal CASSINI experiment of 40 days with 
Sn{f) given by Eqs. ( |4.12 - 4.14 ). The value of the chirp mass A4 ~ 10"^ (for a low 
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frequency cut-off of 10~^ Hz) sets the transition between the systems that are more con- 
venient to be searched during the in-spiral phase and those during the merger phase: for 
M. > 10'' M0, Pmerg rapidly exceeds pinsp by a factor ~ 10 or more. Assuming a low fre- 
quency cut-off around 10~^ Hz, as done in past experiments, the transition occurs at a lower 
value of ~ 5 X 10^ and the drops of (pmsp/Pmerg) is much steeper. As a consequence, 
for Ai ^ 10^ Mq the neglect of the merger signal in searching templates is not serious; for 
Ai ^ 10^ M0 coherent search techniques that include merger waveforms would rapidly in- 
crease the instrument range of sight, by more than one order of magnitude. However, at 
present, the radiation produced during the final burst does not seem to be the most promis- 
ing signal to be searched for in Doppler experiments, as astrophysical considerations suggest 
that its amplitude is not large enough to dominate the noise and matched filtering techniques 
would fail to enhance the SNR due to the poor prior knowledge of the waveform. Only coin- 
cidence experiments can offer a chance of carrying out such searches; this work is currently 
in progress - under the assumption of very simple shapes of the waves - for the analysis 
of data recorded in coincidence during the tracking of ULYSSES, MARS-OBSERVER and 



GALILEO PI 



A. Sensitivity function 

In this Section we discuss the dependence of the sensitivity function T{C,b,C,e,9;aj) on 
the main parameters that characterize the in-spiral signal and the instrument, in particular 
the conditions under which T attains a value close to the ideal one or, equivalently, the 
possible causes of degradation. 

To illustrate the first reason of degradation - the unfavorable angle 6' - as well as the 
behavior of T depending on the frequency range swept by the signal during the observation 
time, we show in Fig. ^ the sensitivity function per unit frequency T(/) for CASSINI's 
nominal noise (see Fig. ^and Table |T|). The sensitivity peak lies in the interval 0.1 ^ ^ ^ 1, 
and broadens and diminishes as cos 6' increases. For ^ <^ 1, T(/) is a monotonic, decreasing 
function of cos 6*, while as ^ exceeds unity, oscillations of increasing frequency occur. The 
overall effect of the angle 6* on T is given by the function T(^e = 0, ^f, = 00, 9; ai), plotted 
in Fig. H, where we also keep ai as free parameter in order to show the effect of the 
redness of the noise spectral density at low frequencies (that is / ^ 1/^); for «! ~ — 0.9, 
T(,^e = 0,^6 = oo,9;ai) has a maximum at the expected value 9 = 7r/2. This is not true 
anymore for redder spectra, and the value of 9 at which T reaches the maximum depends 
on ai; for ai = —2 it corresponds to cos 6' ~ 0.68. 

Concerning the second reason of degradation of T - the limited bandwidth of the signal 
- we show first (Fig. |^) the effect of fb and fe- The deterioration can be strong: T can 
significantly deviate (by one order of magnitude or more) from the ideal value. When both 
and are ^ 1, T is a smooth function of fb and fe, but as soon as they enter the region 
/ ^ 1/T, there are oscillations. We have used the 'nominal' CASSINI form of the noise 
spectrum (Fig. ^ and we keep it fix from now on. 

It is interesting to investigate T in relation to relevant astrophysical parameters. We 
concentrate on four fiducial "best" sources defined by choosing the initial frequency fb and 
assuming that the coalescence occurs at the end of the record, so that = fs and fe = /isco- 
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Each of them is a binary with equal masses (77 = 1/4) , which determines the chirp mass 



M=Mc 



8/3' 



3/5 



(4.27) 



in terms of the critical mass for a nominal run of 40 days, see Eq. ( [4.20[ ). This sets also the 
final frequency /e = /isco and the number of wave cycles recorded at the detector. Table pi 
summaries the main properties of the fiducial sources. In all four cases fe'^ fb, so that we 
have a wide band search and the source's chirp mass is only slightly smaller than the critical 
value. 

Keeping fi, and M. fixed, we first study the effect of Ti on T. As Ti decreases, the signal 
becomes of class III and we loose a progressively greater part of its end. Insisting that the 
chirp mass remains the same, the final frequency fe decreases according to Eq. ( |2.15D . This 
progressive restriction of the search bandwidth is responsible for the great deterioration in 
T, and therefore in SNR, shown in Fig. 10. Since we have assumed that there is no signal 
beyond the coalescence frequency /isco? the sensitivity function remains constant for Ti > 40 
days. 

In Fig. |Tl] we show how the sensitivity function is affected by the round-trip light-time 
T, for the same fiducial sources (and with a run of 40 days) . The variation of T deplaces 
the three-pulse response function - Eq. ( |3.9| ) and Fig. - in relation to the search band 
ifb, fe)- The deterioration with increasing T is due to the fact becomes smaller than unity 
and one looses signal from the high frequency side. 



Finally, in Fig. we investigate the effect of fb (that is the epoch of the binary life at 
which the instrument is "turned on") on T; for each value of /;,, we have in fact computed 
the corresponding critical chirp mass ( |4.20| ) for a 40-day experiment and then evaluated the 
sensitivity function for three different source's orientation in the sky; essentially, for each fj, 
we tune the observation on the mass that produces the highest SNR. The largest values of 
T are achieved for 10"^ Hz < < lO^^ Hz, corresponding to 10^ Mq < Mc ^ 10^ Mq. 



B. Sensitivity of past and future Doppler experiment 

In this Section we analyze the sensitivity of Doppler experiments - in particular CASSINI 
- to search for in-spiral signals; we estimate the range of sensitivity and the achievable signal- 
to- noise ratio as function of the chirp mass Ai and the initial frequency (as well as rj). 
We consider in detail two targets, the Virgo Cluster {D ^ 17Mpc) and the Galactic Centre 
(D ~ 8 kpc); the results are presented so that they can be scaled to an arbitrary distance 
and location in the sky. 

CASSINI will perform three experiments, each lasting 40 days at the epoch of three 
solar oppositions, i.e. when the Sun, the Earth and the probe are approximately aligned, 
in this order, so as to minimize the noise due to interplanetary plasma (see Table | and 



Fig. |T3D [0. The main tracking will be done in X band from the operational antennas 
of NASA's Deep Space Network (DSN), and in band from the new experimental DSN 
station DSS25 at Goldstone (California). Another advanced station to be build in Sardinia 
(the Sardinia Radio Telescope) has been recently funded and it is hoped that it will track 
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CASSINI in Ka-band (up and down) and X-band (down). The radio system of CASSINI's 
mission has a specification for the overall Allan deviation in band of 3 x 10~^^, for 
integration times between 1000 and 10,000 sec, one-to-two orders of magnitude better than 
previous experiments. In addition, the sensitivity will be improved, with respect to past 
missions, by the longer time of observation and, possibly, by the combination of the three 
records. Fig. |13| shows the round-trip light-time and the value of cos 9 for the Virgo Cluster 
and the Galactic Centre during the cruise of the space-probe to Saturn. 

We begin by investigating the expected SNR produced by in-spiraling signals of possible 
sources located in the Virgo Cluster in the plane {A4,fb)- Since it turns out that the 
three planned experiments have comparable sensitivities, we give explicit results only for 
the last one, in 2004. Fig. |14| shows for the Virgo Cluster three different levels of SNR: 
p = 1 , 5 and 10. It confirms the theoretical suggestion that the detectable events of class 
III lie in a narrow strip below the critical value of the chirp mass. All events above it are 
of class IV and include the merger. This plot shows how the SNR changes as function of 
both A4 and /;,: at p = 1 and for black holes of comparable masses, the smallest detectable 
chirp mass is ^ 5 x 10^ Mq. Lower masses are excluded because the systems are too far 
from coalescence and masses larger than ~ 10^ Mq are not accessible because their emission 
frequencies are too low. Note the drastic curtailment if a low frequency cut off at 1/T ~ 10~^ 
Hz is assumed, as it was done in the past with ULYSSES; the strongest signals fall below 
ffyT = 1. The SNR is very much reduced when a small mass ratio is assumed (lower panel); 
an important and unknown factor in assessing the detectability of gravitational radiation 
from MBHB's is the statistics of this ratio. 

One wonders if a source in the Virgo cluster detected in the last opposition is traceable 
also in the second or the first (albeit at a smaller SNR), thus adding important phase 
information. We checked that for Ti = 40 days this is not the case: detectable sources are 
too near coalescence, in relation to the (roughly) one year separation between different runs. 
We have then tried doubling the observation time to Ti = 80 days ('extended' experiments) 
and interestingly found that for some sources the three extended data sets can be combined: 
the strongest sources detectable in one opposition are also visible in the previous one. A 
single, extended experiment does not show drastic improvements over the nominal record; 
however, the accessible portion of the plane {J^,fb) becomes slightly wider. Of course, 
a drastic improvement in sensitivity and detection probability would be expected for a 
continuous run of a few years. 

How far is the farthest detectable system? From Fig. we see that the highest SNR, 
about 28, produced by a source in the Virgo cluster is attained for a system of comparable 
masses ~ 5 x 10^ Mq; at p = 1 this corresponds to a distance ~ 600 Mpc; if M2/M1 = 10~^, 
the ken is reduced by a factor ~ 5. 

In order to have an idea of the change in SNR due to different angles 6, in Fig. |16| we 
have studied this dependence for the third experiment in three different cases: Ai = 10^ , 10^ 
and 10^ M0 (/ft = /s, = /isco and t] = 1/4). Note that, for M = 10^ Mq, cos9 = is a 
local minimum, rather than the expected maximum: this is due to the fact that signals fall 
in one of the high frequency minima of the three-pulse response, see Figs. § and 0. Figs. |l^, 
p!5| and |T^ therefore allow to have an overall picture of the all-sky sensitivity of CASSINI, 
see Eqs. (|A6]) and (|4T7| ). 

We turn now to galactic observations and consider a binary system in the centre of Our 
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Galaxy composed of a primary black hole of mass Mi = 2 x 10^ Mq, whose existence is 
strongly inferred from observations carried out over the past 25 years and a secondary, 
smaller object of mass M2. In Fig. |1^ we show for the third CASSINI opposition the 
SNR contour levels p = 1,5, 100 in the plane (M2, Of course, there is no signal above 
/isco ^ 2 X 10^3 (Ml/2 X 10^ Mq)-^ Rz, dr. Eq. (gl|), essentially independent of M2. 
CASSINI proofs to be a fairly sensitive instruments for such signals; in fact, each of the 
three experiments is sensitive to secondary objects of mass down to ~ 10 — 50 Mq and 
for M2 ^ 100 Mq the expected SNR can be ^ 10, see Fig. |1^. Remarkably, most of the 
sources visible during one experiment are detectable also in the others; in fact the time to 
coalescence is ~ 2995 {fb/5 x 10"^ Hz)-^/^ ^iqq Mq/p) (2 x 10^ Mq/MY'^ days (cfr. also Fig. 
191) . If a candidate signal in spotted in one of the data sets, it is therefore possible to chase 



it also in the other oppositions, which provides a robust method of accepting or discarding 
candidate events by coherently tracking the gravitational phase over a period of some years. 
This is a direct consequence of the fact that many signals would be recorded as class II 
and III (therefore far from the final coalescence) as is shown in Fig. |1^ by the dotted and 
dashed line in the plane {M2, fb)] this is in contrast with what appends for observations 
in the Virgo cluster, where most of the signals show up as class IV. The good sensitivity 
for nearby sources allows also to search for possible binary systems in other galaxies of the 
Local Group. 

It is interesting to compare CASSINI's sensitivity with past Doppler experiments (Table 
I). The best Doppler data available today come from observations carried out in 1991, 1992 
and 1993 with the probes ULYSSES, GALILEO and MARS-OBSERVER. In 1993 all three 
spacecraft have been simultaneously tracked for about two weeks, with a triple coincidence 
experiment whose data are being analyzed; this experiment will be useful to search for 
signals coming from the merger phase . 



We concentrate on ULYSSES' 1992 experiment |T^,|T5|. After pre-processing, only 14 



days of data have been retained, centered around the sun opposition. Their Allan deviation 
(at 1000 sec) is variable, with only its best value near the specification, 3 x 10~^^. Data 
outside the bandwidth ~ 2.3 x 10~^ — 5 x 10~^ Hz have been discarded, see Fig. ^: at high 
frequency they are corrupted by the anomalous rotation of the spacecraft; at low frequency, 
the subtraction of orbital and media contributions was carried out with a fit, which made 
the data unreliable below / ^ 1/T. In the intermediate band the noise spectrum is well 
approximated by a power- low oc /~°'^ (see Table 0). The low frequency cut off limits the 
largest observable mass to 

< 9 X 10^ {4r]f^ Mq , (4.28) 

and the value of the critical chirp mass ( [4. 201 ) for the 14-days ULYSSES second opposition 
experiment is 




^ 1.3 X 10^ -^zrrr — 1— ^0 ? (4-29) 



Eqs. ( [4.28| ) and ( [4.29| ) clearly indicate that the most interesting region in the {Ai, fb) plane 
is not accessible (cfr. Fig. |1^) and therefore the largest attainable distance is strongly 
reduced. Moreover, the noise spectral density, as defined by Eqs. ( |4.12| )-( |4.15| ) with the 
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parameters reported in Table H (see Fig. is about two orders of magnitude larger than 
for CASSINI. According to Eq. (^4.101) , the ideal ken (we recall that it is defined for p^^ = 1) 
holds 

^■^"H i-3x10hJ ( 2.5 X 10-^^ Hz- ) (iSsO^) 

Therefore, in conventional astrophysical scenarios, ULYSSES' experiments are relevant only 
for the Galactic Centre, whose angular position with respect to the spacecraft was 6 ~ 109". 
It is straightforward to check that at p = 5 secondary black holes with M2 ^ 10^ Mq would 
be detectable. 



V. CONCLUSIONS 

We have analyzed the sensitivity of Doppler detectors to gravitational waves emitted by 
coalescing binaries in their in-spiral phase. For CASSINI's experiments sources in the Virgo 
Cluster would be observable at a signal-to-noise ratio up to ~ 30. Furthermore, binary 
systems with Mi ~ M2 ~ 5 x 10^ Mq at distance ~ 600 Mpc are within the range of the 
instrument. If in the centre of Our Galaxy there is a massive black hole of 2 x 10^ Mq, 
orbiting black holes with mass ^ 50 Mq would be easily detectable; for larger masses the 
signal-to-noise ratio can be considerably higher and therefore one can reach other galaxies 
within the Local Group. These are considerable improvements in sensitivity and range of 
accessible masses with respect to past Doppler experiments. 

This analysis has been made in relation to astrophysically reasonable sources, but without 
any concern for the probability of such events and their detection. Considering our ignorance 
about sources of gravitational waves, in particular for events in galactic nuclei, this point of 
view appears reasonable, especially for non dedicated experiments like CASSINI's; however, 
it is healthy and sobering to assess what one can predict on the basis of very simple and 
conventional astrophysical models. The estimation of the event rate is a much debated 



problem in low frequency gravitational wave experiments (see |T^,|32|J3^,|36|J3^ and reference 
therein), with sensitivity not good enough to detect signals from known stellar mass galactic 
binaries. We briefly consider here two models: (a) accretion of small black holes onto a 
massive one in a galactic core and (b) coalescence of two MBHB's of comparable mass, 
possibly resulting from the mergers of two galaxies. In the context of Doppler experiments, 
they are appropriate to small and large distances, respectively. 

Let consider now case (a). Suppose that a central black hole M in a typical galaxy 
(i.e., in Our Galactic Centre) has captured = M/M2 masses M2 at the uniform rate 
R = H0M/M2 during the whole life of the Universe H^^ ~ 10^°/iiooyr, where /iioo is the 
Hubble constant normalized to 100 kmsec"^ Mpc~^; each of these captures gives rise to a 
chirp. In this ideal case, in which the spread in mass is neglected, one could consider the 
search as follows. For a given M2 (and fixed Mi, which determine uniquely the chirp mass) 
there is a characteristic initial frequency Jb, given by Eq. ( p.l8|) , for which we have a wide 
band detection and coalescence at the end of the record; this corresponds to the situation in 
which the SNR is maximum. As diminishes, the fractional bandwidth A - see Eq. ([4.21 
- decreases according to 
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and so does the SNR, proportionally to A: 



A)-8/3] 



3/8 



P 



(5.1) 



(5.2) 



where f ^ fb ^ fe- Since a small value of A increases the useful time for detection Td - that 
is the total time for which the source is "visible" above the detection threshold po ^ detection 
will generally occur at a SNR less than the maximum (attained for fb = Jb and fe = /isco) 
and with a start frequency < fs- po determines the threshold Ad for the fractional 
bandwidth and Td, via Eq. ( |4.22|) . Indeed, the success probability is P = RTd and we have 
computed it using the results reported in Fig. |T^. Its behavior as a function of the mass 
M2, which enters through both R and Td, is given in Fig. |19]for the thresholds pd = 1? 5, 
10 and 15. This allows us to determine the most favorable value of M2, about 100 Mq, and 
the corresponding (optimistic!) success probability, in each experiment of Ti = 40 days, is 
~ 10-^ 

We turn now to case (b). To go to very large distances and to reach the largest number 
of sources, we must aim at large masses and use low frequencies, in the band fT^l. In 
this frequency range, the noise spectral density is well described by a single power low 



Snif) ^ So (jr 



(5.3) 



The SNR produced by a source that sweeps the bandwidth A, where = (1 + ^) ^b, see Eq. 
( [4.21| ), can be analytically computed from Eq. (|4.16| ) assuming the simple low- frequency 
behavior of the Doppler response filter ( p.lO|) : 
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is a numerical coefficients that depends on the angles t , and 9. The largest attainable 
distance, in observations of a binary of mass Ai is therefore 
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and the corresponding number of detectable events 
N ^RTi {DHof 

- PT- u-3 ^3/2 M^^"^ p3ai/2 .l-3«i/2 
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Of course, the range of the instrument depends also on the mass of the source and one can 
express as a function of /;,, Ti and A, using Eq. ( |4.19|) ; Eq. (|5.6|) therefore becomes: 
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where k' is a numerical constant depending on k. We see that N is proportional to ^ ' 

and increases with decreasing frequencies if a\ < —2. The nominal case ai = —2 (Table 
II) marks the transition between the two alternatives, whether it is useful or not to go 
to lower frequencies. If we assume (naively and optimistically) that each galaxy in the 
Universe contains a massive black hole and has undergone merging events leading to 
the formation of a binary with time to coalescence shorter than the Hubble time, we can 
estimate the event rate as R = N^Ng/ Hq, where Ng is the total number of galaxies. As we 
have seen in the previous section, the ken of CASSINI is D ~ 600 Mpc and therefore the 
probability of detection is 

Of course, larger distances could be attained if we knew the waveform produced during the 
final collapse and the merger could be included in the search. In the present situation of 
ignorance, however, a good level of confidence in the detection of merger signals requires 
coincidence experiments. 

We conclude that, in conventional astrophysical scenarios, the events we are looking 
for do not have a large probability to be detected, the extragalactic case (b) being more 
favorable. Note, however, that for low chirp mass sources A4 ~ 100 M©, the total SNR 
could be increased by combining two or three records. 

Present Doppler experiments are relatively cheap, but use non dedicated spacecraft and 
are limited in instrumentation, orbits, observation time and by several other constraints. 
However, despite their low sensitivity, they are the only one that at present can access the 
coalescence of binaries involving massive objects, and therefore quite valuable. They might 
discover gravitational waves; if they do not, they will provide interesting astrophysical limits 
on massive sources, in particular in Our Galactic Centre, the Local Group and the Virgo 
Cluster. Since real data are already available in large amounts and the search techniques 
are very similar to those implemented for Earth-based (and, in future, space-borne) laser 
interferometers, Doppler experiments can also be used to test filtering and processing tech- 
niques. 
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TABLES 



TABLE I. Summary of Doppler experiments to search for gravitational waves. Ti is the 
effective length of data available; ay the average effective Allan variance of the data and T the 
average round-trip light-time of the radio link during the data acquisition runs. The frequency 
bands of the up-link (S: 2.1GHz, X: 7.2 GHz, K„: 34 GHz) and down-link (S: 2.3 GHz, X: 8.4 GHz, 
Ka'. 32 GHz) carriers are shown in the two final columns. For the four short experiments carried 
OTit. witli PIONEER 10 and 11 and VOYAGER 1 wx^ [21] and ix^'tnxnicc^s ihcivm. 



Space-probe 


Ti/days 


y 




T/sec 


up-link 


down- link 


ULYSSES 3-Dec-90/4-Jan-91 


3.5 


3 X 10" 


14 


600 


s 


S,X 


ULYSSES 20-Feb-92/18-March-92 


14 


7x 10- 


14 


4430 


s 


s,x 


ULYSSES Mar-93/Apr-93 


19 


1.4 X 10 


-13 


4100 


s 


s,x 


GALILEO Mar-93/Apr-93 


19 


2.3 X 10 


-13 


945 


s 


s 


MARS-OBSERVER Mar-93/Apr-93 


19 


5 X 10- 


14 


1150 


X 


X 


CASSINI 26-Nov-01/5-Jan-02 


40 


3 X 10- 


15 


5783 


X, Ka 


X, Ka 


CASSINI 7-Dec-02/16-Jan-03 


40 


3 X 10- 


15 


7036 


X, Ka 


X, Ka 


CASSINI 15-Dec-03/24-Jan-04 


40 


3 X 10- 


15 


7859 


X, Ka 


X, Ka 
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TABLE II. Characteristic parameters of the ULYSSES and CASSINI noise spectral den- 
sity. The table shows the reference frequency 1/T; the reference noise spectral density level 
So = Sn{l/T); the transition frequencies /i2,/23; the exponents aj and the frequencies fj 
(j = 1,2,3) of each regime, see Eqs. ( 4.12| ) - ( |4.14 ). When, as in our case, r/12 = 1, r/23 » 1 
and as > 02 then T/i = 2V"i, T/2 = 2V°2 and T/3 = r/23 2i/"3 (r^s)-"^/"^. 



parameters 


ULYSSES 


CASSINI 


(i/r)/Hz 


2 X ID 


10 


5o/Hz-i 


2.5 X 10-24 


9 X 10-26 


/12/Hz 


2 X 10-4 


10-4 


/23/Hz 


5 X 10-2 


10-1 


/i/Hz 


2 X 10-4 


7.07 X 10-5 


/2/Hz 


8 X 10-4 


2.5 X 10-5 


/3/Hz 


5 X 10-2 


7.95 X 10-1 


Ql 


-10 


-2 


02 


-1/2 


-1/2 


"3 


10 


2 



TABLE III. Main properties of the four fiducial sources used in the text. The assignment 
of the initial frequency fb, in the case of equal masses and coalescence at the end of the record, 
determines the final frequency fe (for a 40 days run), the chirp mass M and the number of wave 
cycles M recorded at the instrument while the signal sweeps the frequency interval {fb, fe)- 



fb/H^ ^/Hz M/Mq 

10-5 5.3 X 10-5 3.6 X lO'^ 52 

5 X 10-5 6.9 X 10-4 2.9 x 10^ 273 

10-4 2.1 X 10-3 9.2 X 105 550 

5 X 10-4 2.7 X 10-2 7.1 X 104 2761 
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FIGURES 




FIG. 1. Classes of coalescing binary signals for one of the three CASSINI's experiments 
(Ti = 40 days). The upper panel shows the partition of the plane {M., fb) in periodic (I), linear (II) 
and non-linear region. The dotted line divides non-linear chirps in signals of class III and IV. The 
dashed lines indicate the boundary between signals of class IV and bursts - i.e. /isco as a function 
of Ai - for M2/M1 = 0.01 (lower line) and 1 (upper line). The lower panel shows the Newtonian 
times to coalescence r„, Eq. ( |2.23 ), computed at the transition frequency (see upper panel) for 
each class of signals. 
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FIG. 2. To show how the partition of signals into different classes depend upon the ob- 
servation time Ti we consider two binary systems: Mi = M2 = lO'^ Mq (left panels) and 
Ml = 2 X 10^Mq,M2 = 1000 (right panels). The Newtonian time to coalescence r„, Eq. 
( 2.23| ), and the initial frequency Z^, Eq. ( 2.22| ), are used to separate the classes: continuous line 
between I and II, dashed-dotted line between II and III, dotted line between III and IV. The 
frequency scale has been chosen so that /isco corresponds to the top of the panel. 
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FIG. 3. The effect of the Doppler filter rg on in-spiral signals. The plots show a comparison 
between the chirp amplitude A{t), Eq. (|2.5D , (solid line) and the Doppler amplitude for different 
values of the angle 6 (cos 9 = 0: dashed-line; cos 9 = 0.4: doted-line) and of the chirp mass: (a) 
M = 10^ Mq, (b) M = 10^ Mq, (c) M = 10^ Mq, (d) M = 10^ M© (r/ = 1/4 for all systems). The 
chirp amplitudes are normalized to A{t = 0) = 1; the round-trip light-time is T = 4000 sec. Note 
the increasing modulation induced by the Doppler response as the final merger is approached; this 
modulation disappears as M./T increases because the filter degenerates in a multiplicative factor 
(cfr. Eq. (|^ and Figl). 
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FIG. 4. The square modulus of the Doppler filter Eq. ( |3.9|) , as a function of ^ = f T 

for selected values of the angle 6 (solid line: cos 9 = 0; dotted-line: cos 9 = 0.4; dotted-dashed line: 
cos^ = 0.8). The small box zooms the behavior of for ^ 1; note that, only when ^ is 

small, is a monotonic decreasing function of cos^, which is not true for ^ 1. 
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10-2^ r 



frequency / Hz 



FIG. 5. Noise spectral density of CASSINI (solid line) and ULYSSES (dash-line). 5„(/) is 
computed according to Eqs. ( 4.12 ) - ( [4.15| ) and Table |l|. 
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FIG. 6. Comparison of the sensitivity of CASSINI experiments in observations of in-spiral 
and merger signals. For a given source, the plot shows the ratio of the optimal SNR achievable 
in observations of the in-spiral signal, pinspi with respect to that produced during the merger 
phase, /Omerg) as a function of the chirp mass and for the mass ratios 1 and 0.01, label (a) and (b) 
respectively. The solid lines refer to the low frequency cut-off 10~^ Hz, whereas the dotted lines 
to 10~^ Hz. The small box shows, for the same cases, the relative overlap J, Eq. ( [4.26| ), of the 
in-spiral and merger signal with respect to the instrument response. The time of observation is 
assumed to be Ti = 40 days and the noise spectral density Sn{f) is computed according to Eqs. 
(P^ ) - (p^ and Table 0. 
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FIG. 7. The spectrum T(/) of the sensitivity function, see Eq. ( 4.17| ). We have used 
CASSINI's noise parameters, Table ^ and Fig. |5| (sohd hne: cos 9 = 0; dotted hne: cos 9 = 0.4; 
dotted-dashed hne: cos 9 = 0.8). 
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FIG. 8. The deterioration of the sensitivity function T, and therefore of the SNR, due to the 
angle 9 between the source and the probe and the redness of the spectrum at low frequencies, given 
by ai (cfr. Eqs. (4.12) - ( [4.15 ) and Table ||), as described by the function T(,^e = 0, = oo, 6\ ai). 
In the upper panel we show, as a function of ai, the maximum value of T(^e = 0,^;, = oo,0;ai) 
(triangles) and the angle ^max (bullets) at which it is attained. The 'normal' value 0max = 71/2 oc- 
curs above ai ~ —0.9. The lower panel shows the behavior of T(^e = 0, = oo, 6] ai) as a function 
of cosO for the case ai = (dashed line), ai = —1 (dotted-dashed line), ai = —2 (continuous line) 
and ai = —3 (dotted line); the values of a2 and as are —1/2 and +2, respectively, corresponding 
to those of the CASSINI experiments. Note that, for qi < —1, the smallest deterioration does not 
occur at the obvious value 9 = 7r/2; in particular, for ai = —2 it occurs at cos9 ~ 0.68, so that a 
small signal enhancement with respect to the ideal case is possible. 
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FIG. 9. The behavior of the sensitivity function T as a function of the final frequency /e (in 
units of the round-trip hght-time T) for selected values of the initial frequency T = ^f, = 0.1 , 1 , 10 
(see labels). CASSINI's noise parameters have been used (Table |l| and Fig. |). To the right of the 
arrows we have broad band searches with fe — fb > fb (solid line: cos ^ = 0; dotted line: cos 6 = 0.4; 
dotted-dashed line: cos^ = 0.8). 
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FIG. 10. The behavior of the sensitivity function T as a function of Ti, for the four 
fiducial sources described in the text and in Table |l| (panels (a) to (d), respectively), with 
fb = 10~^ , 5 X 10^^ , 10^^ , 5 X 10^^ Hz, equal masses, coalescence at the end of the record and 
chirp mass almost equal to its the critical value for a 40 days run. The deterioration as Ti de- 
creases reflects the loss of signal-to-noise ratio due to the narrowing of the band-width (solid line: 
cos 9 = 0; dotted line: cos 6 = 0.4; dotted-dashed line: cos 9 = 0.8). Notice the different scale in 
panel (a). 
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FIG. 11. The behavior of the sensitivity function T as a function of the round-trip hght time 
for same four fiducial sources of the previous Figure (sohd Une: cos 9 = 0; dotted hne: cos 6 = 0.4; 
dotted-dashed hne: cos 9 = 0.8). Notice the different scale in panel (a). 
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FIG. 12. The behavior of the sensitivity function as a function of the initial frequency ft, for 
wide band searches of class III, with ft = Ib, fe = /isco and Ti = 40 days. The corresponding 
value of the chirp mass - near to the critical value ( [4. 201 ) - is given by the filled circles (solid line: 
cos^ = 0; dotted line: cos 9 = 0.4; dotted-dashed line: cos 9 = 0.8). 
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FIG. 13. CASSINI's cruise to Saturn. The solid line gives the round-trip light-time T as 
function of time and we show also the variation of cos 6 for the Virgo Cluster (dashed-dotted line) 
and the Galactic Centre (dotted line). The three gravitational wave experiments take place during 
the days 695-735, 1071-1111 and 1443-1483 from 1st January 2000. 
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FIG. 14. Detectable in-spiral signals emitted by binary systems in the Virgo Cluster for the 
third CASSINI's experiment of 40 days. The contour plots of the SNR (in amphtude) for the levels 
1, 5 and 10 are shown. The top and bottom panels refer to the mass ratios M2/M1 = 1 and 0.01, 
respectively. In the latter case no signals are observable at p > 10. The dashed diagonal line gives 
the critical chirp mass see Eq. ( |4.2C| ). The dotted line marks the transition between signals 
of class III and IV. Notice that for Mi = M2 (upper panel) the two lines almost coincide. 
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FIG. 15. Signal-to-noise ratio for in-spiral signals emitted by binaries in the Virgo cluster 
(D = 17 Mpc) with coalescence at the end of the record - that is fb = fs and /e = /isco - as a 
function of the chirp mass. The solid line refers to the case Mi = M2, while the dashed line to 
M2/M1 = 0.01. We have considered CASSINI third opposition and Ti = 40 days. 
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FIG. 16. The signal-to-noise ratio as a function of cos Q during CASSINI third opposition, 
T\ = 40 days, for in-spiral signals emitted by binary systems located at the same distance of the 
Virgo cluster {D = 17Mpc) with /{, = and /g = /iscoj that is final coalescence at the end of the 
record (squares: M = 10^ M©; triangles: Mq = l^'^ Mq- bullets: Mq = lO^M©; for ah = 1/4). 
The decrease for | cos^| — > 1 is of course due to the approach to the case of propagation along the 
instrument's beam; the saddle in correspondence of | cos 0| ~ for Ai = 10® Mq is due to the fact 
that the frequency band swept by the signal encounters a minimum of the response function at 
fT ~ 1, which vanishes for cos^ = 0, see Figs. Q and |^. 
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FIG. 17. Detectability of radiation emitted by the in-spiral of a secondary black hole of mass 
M2 orbiting a primary object with mass Mi = 2 x 10^ M©, in the Galactic Centre (D ~ 8 Kpc) 
during CASSINI's third experiment (Ti = 40 days). The contour plots of SNR = 1 ,5 and 100 in 
the plane {M2, fh) arc shown. The dashed-dotted line divides the plane (M2, ft) in signal of class 
II (below) and III (above), while above the dotted line the signals are registered as of class IV. 
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FIG. 18. Expected optimal signal-to-noise ratio produced in CASSINI records by in-spiraling 
binaries in the Galactic Centre. The plot shows the maximum SNR - i.e. computed for = Jb 
and coalescence at the end of the record (Ti = 40 days), therefore fe = /isco - produced by binary 
systems with Mi = 2 x 10^ Mq and a secondary black hole of mass M2 < Mi as a function of M2, 
during CASSINI's third opposition. 
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FIG. 19. Time of 'visibility' To and probability of success for the detection of radiation 
emitted by the in-spiral of a secondary black hole into a primary of mass Mi = 2 x 10^ Mq in Our 
Galactic Centre during CASSINI third opposition. The curves refer to four different thresholds of 
detection (pD = 1: solid line; po = 5: dotted line; po = 10: dashed line; pD = 15: dotted-dashed 
line) . 
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